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The study of the topology of energy bands in sohd has always been interesting 
and fruitful. Historically, Thouless et al proposed the TKNN number or Chern 
number of the energy band to explain the quantization of Hall conductance in 
the integer quantum Hall effect. Recently, Z2 topological insulators have been 
intensively studied. There are gapless edge states in these materials which are 
protect by bulk band topology. The topology of the bands is characterized by 
some Z2 invariants. However, Z2 invariants are crude and strongly dependent 
on the symmetry. Here we give an unified picture of the relationship of the 
edge states and the geometry of the energy bands. We show the band geome- 
try determines not only the topological but also the geometrical properties of 
the edge states. Our picture is applicable even when the symmetries are broken. 



To explain the quantization of Hall conductance in the integer quantum Hall effect 
Thouless et al calculated the Hall conductance of the two dimensional electron gas in a 
uniform magnetic field and period potential using the Kubo formula[l|. It's shown the 
quantization of Hall conductance is a bulk property and directly related to a the TKNN 
number or Chern number of the energy bands. However, the quantum Hall effect in a 
finite system with edge depends on the existence of the gapless edge state. Then Hatsugai 
established correspondence between Chern number and gapless edge states j2|. In 1988 
Haldane proposed a tight-binding model which support a occupied energy band with 
nonzero Chern number in ground state [3]. The system can exhibit quantum Hall effect in 
the absence of macro magnetic field. We refer such materials as Chern insulators. The 
Chern insulators only exist when time reversal symmetry is broken. Some of the early 
proposed time-reversal symmetric topological insulators can be viewed as two copy of 
Chern insulators. In these insulators there are one component of the spin or pseudospin 
is conservedj4-7|. So the occupied bands of the ground states can be easily divide to two 
bands with opposite nonzero Chern numbers according to the value the spin or pseudospin 
takes and there are clear energy band and edge state correspondence just like in the Chern 
insulators. In these cases the physics pictures are quite simple: the two time-reversal related 
bands pump the spins or pseudospins to opposite directions when electric field present 
and the quantum spin Hall effect can be easily interpreted in this way. However when 
there are no conserved spin or pseudospin as in js] when the Rashba term are nonzero, the 
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two occupied bands don't naturally form a time-reversal related pair. So it's claimed the 
Chern numbers lost meaning and the Z2 topological invariant is proposed to characterize 
the statejol. The classification is mathematically involved and make the physical 
picture complicated. It's even claimed the isolate Z2 insulator will return to its origin 
state after two pumps 1^ . We will show the Chern number can still be used to classify 
the two dimensional topological insulators and the Z2 topological insulator can be viewed 



as two Chern insulators put together 



are proposed and observed shortly after 



The three dimensional topological insulators 

n 

13H15|. It's classified by four Z2 invariants. The 



strong topological insulators have odd number of Dirac points at the surface. The existence 
of these Z2 invariants depends on the time-reversal symmetry. With the similar idea a 
topological crystalline insulators are proposed with C4 or Cg symmetry and the Z2 invariant 
characterize this insulator depends on these symmetries 15| . 



All those topological invariants are somewhat crude, that is they only determine the 
topology of the edge state. Historically the word "topology" was first used to distinguish 
"qualitative geometry from the ordinary geometry in which quantitative relations chiefiy are 
treated" (if)!] . The quantitative geometry certainly contains more details than the "qualitative 
geometry". In this article we show how the geometrical property of the energy bands is 
studied and how it's related to the edge state property. The first problem is to distinguish 
individual energy band. In most of insulators there are several occupied energy bands, e.g. 
because of time-reversal symmetry ground state of two dimensional Z2 insulator have even 
number of occupied bands. These bands form a trivial vector bundle. The topological 
non-triviality come from the time-reversal symmetry. We can always find global continuous 
bands that span a trivial bundle. However, if we insist on using two groups of time-reversal 
related bands to span the bundle, the two groups of bands won't be global continuous and 
each has an odd Chern number [HI]. We propose a natural way to define the two group of 
time-reversal related bands. In this way we can establish the energy band and edge state 
correspondence just like in Chern insulators. Let ai and a2 denote the primitive vectors 
of the two dimension lattice, ai parallel to the edge. The wave vectors of the occupied 
bands are functions of ni = k ■ ai and 71,2 = k ■ a2. For simplicity we assume two bands 
are occupied, all our results can be easily generated to more bands. Two global continuous 
orthonormal wave vectors ]'Ui(rii, 77-2)) and ]ti2('^i) ^2)) can be defined to span the occupied 
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space. When ni is fixed the problem become one-dimensional. To make a natural decoupling 
of the occupied bands in one dimension we demand the in each of the decoupled bands the 
wave vectors are smoothly defined, that is they can be connected by parallel transport along 
72-2. Parallel transport is a geometrical property of the occupied bands and the geometrical 
structure of the occupied bands are determined by the non-Abelian Berry connection. From 
real space point of view the natural decoupling should make the one-dimensional Wannier 



functions as localized as possible. In fact the two approaches give the same result 



19|. A 



parallel transport of wave vectors \ui){i = 1, 2) at n2 = n2o along n2 can be expressed as 



{\vi{n2)),\v2{n2))) = {\ui{n2)),\u2{n2))) Pexp [i A2{n2)dn2] (1) 

'n20 



"2 



where 



{U2\£;\ui) (W2|J;|W2) 

is the non-Abelian Berry field, it's a 2 x 2 matrix 17|. For notation convenience we drop 
the ni temporarily. \vi{n2)), \ v2{n2)) are smooth along n2. However there is a mismatch 
Ug{n2o) = Pexp (^i /^^^J'^^'' ^42(^2)^^2^ between \vi{n2o + 2%)) and the \vi{n2o) = \ui{n2o). 
To cancel this mismatch we diagonalize U{n2o) by a U{2) transformation. Then we 
perform the same U{2) transformation at \ui) for all 77.2 and in this way we establish a new 
frame. It's easy to prove in the new frame there are only two U{1) mismatches between 
\vi{n2o + 27t)) and the \vi{n2o)) = \ui{n2o)), which are the eigenvalues of the Ug{n2o)- In 
the new frame \vi{n2)) can be connected by parallel transport except an U{1) number 
and the occupied bands are naturally decoupled to two. The eigenvalues exp{i2n^i) of 
Ug{n2o) are independent of 7720 and if they are nondegenerate the decoupling will be unique. 
For all 77i we choose 7720 = and make the decoupling. Here we only consider the case 
eigenvalues of Ug{n2o) are degenerate for isolate 77i and the decoupling at the degenerate 
points can be determined by continuous extension. In this way the two occupied bands of 
the two-dimensional system are naturally decoupled and the bands are smooth along the 
72,2 direction. All our discussions are based on such a decoupling which is certainly depend 
on the particular edge we study. So we don't study the general geometry of the bands but 
the geometry of bands for a given edge. defined modulo 1, give real space position of 
the center of the Wannier functions. $j exhibit the geometry property of the bands. The 
change of $j in one period of 7ii(27r) is the Chern number of the band. In the Z2 insulator 
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because of the time-reversal symmetry the decoupled bands are related by time-reversal 
operator and = $2(— ''^i)- The two bands have opposite Chern numbers and Z2 is 



the parity of the Chern number. In [18| with a simplified model of edge topology of the 
edge state is identified with that of centers of the Wannier functions So we have the 

energy band edge state correspondence: each of the two dimensional naturally decoupled 
bands with nonzero Chern number determines a branch of gapless edge state. In two 
dimension insulator the study of topology of energy bands is reduced to study the topology 
of the lines formed by $iS. Same idea can be applied in three dimension. In three dimension 
the topology of the bands are determined by the topology of surfaces of $j(ni,n2), where 
Tix = k • ai,n2 = k ■ a2 and ai, a2 are the primitive vectors at the two dimensional edge. 
However, in three dimension, because of the complexity of the topology of the surface, e.g. 
existence of the "Dirac" type of points, the energy bands can't always be decoupled. The 
above discussions can be easily applied to the time-reversal symmetry broken materials 
with several occupied bands. In this work we study several types of topological insulators 
with and without symmetry protection by direct comparing the geometry of the $j and the 
spectrums of edge states. We conclude that the geometry of the naturally decoupled bands 
determine not only the topology but also the geometry of the edge state. So the edge state 
of those topological insulators actually protect by the geometry the occupied bands. 

Our first example is graphene, the first proposed Z2 topological insulator jsj. In FIG 1 we 
compare the geometry of the $iS and that of the spectrum of edge states. We see edge state 
occurs at where the geometry of the bands is more nontrivial and there are no edge states 
at where $jS are relatively flat. The geometrical structure of the spectrum of the edge state 
and that of <l>j are very similar. When the parameter A^, varies granphene experience a 
topological and trivial insulator phase transition. The $j not only indicate the topological 
change but also accurately predict where the gap of edge state opens. Similar result can 



be obtained from the HgTe quantum well. We use the simple tight-bonding model in jo] 
with the d = 70A. We can see $jS are quite fiat except a small interval around the ni = 
where the the edge state occur. We also calculated the $jS of a Z2 trivial topological 
insulator in |^]. The occupied bands are naturally decoupled to two time- reversal related 
bands with Chern number ±2. Compare the $jS with the edge state spectrum we see the 
gapless edges states are actually protected by the geometry of the bands. In this case 
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the degeneracies at the cross points of the edge state spectrum are very similar to the 



"accident degeneracies" introduced by zone folding when cell are enlarged [20|]. It's claimed 
these degeneracies aren't protect by symmetry so can be easily lift by random potentials 
lower the translational symmetry. However numerical simulation shows they are at lest 
quite strong against time-reversal symmetric random potentials along the n2. Though the 
position of the cross points varies with different random potential the degeneracies aren't 
lifted by rather strong potentials. So whether the geometrical protected edge states in this 
case is trivial is still an open problem. 

In Fig 2 we show some comparisons between the $jS as function of rii and 72-2 which 
are surfaces and the spectrum of edge state of the tree dimensional topological insulator 
in 15| (For convenience the we draw the energy spectrums which are most close to the 
bulk gap, only the part cross the bulk gap are actually are edge state). We see the 
conclusions from two dimensional case are still applicable here: edge states only occur 
when the geometry of $jS are more nontrivial. Moreover, in the cases of three dimensional 
topological insulators our geometrical representation shows great advantage over the 



topological indices in 15|. The geometrical representation predicts not only the existence 
of the edge state Dirac cones but also where the Dirac cone appears. So the topological 
indices only give part of the topological information that can be inferred for the geometrical 
structure of the surfaces of the $iS. Our geometrical representation clearly show the 



picture of three dimensional topological in 2M, |22| fails. Their picture claim there must 



be monopole of Berry field for the Chern number change between different surfaces in 
Brillouin zone which means occupied bands must "collide" with each other at the middle 
of half Brillouin zone. From Fig 2 we can easily observe that the differences of Chern 
numbers of the bands in the two dimensional surface come from the switch partners at 
the Dirac type of points of ^iS. In Fig 3 we show comparisons <l>jS and edge state spec- 



trum in Be2Se3 with the tight-binding model in 23|. The correspondence are almost perfect. 



In 15| the so called topological crystalline insulators are proposed and it's claimed the 
gapless edge states in these materials are protected by the time reversal symmetry and 
the C4 or Cq rotation symmetry. A rather complicate Z2 topological invariant is proposed 
to characterize the topology od the occupied bands and the invariant depend on both the 
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time reversal symmetry and the C4 or Cq rotation symmetry. It's shown in Fig 4 the $iS 
not only indicated the topological property of edge state spectrum but also accurately 
predicted the geometrical structure of the edge state. We also deliberately break the C4 
rotation symmetry by using different parameter t'2 in x and y directions in the proposed 
tight-binding model. In this modified model there is only a C2 rotation symmetry left, the 
discussions of the topological properties of the bands in |l5| are no longer applicable, that 
is, the Z2 invariant can't be defined here. However, we shows the edge states geometry or 
topology aren't changed by the broken C4 symmetry even with a very large difference of 
in X and y directions. In Fig 4 we see the $jS can still predict the edge state geometry when 
the symmetry is broken. So we propose in this case it's more convenient to use the geometry 
of the bands to characterize the insulator than to use some complicated topological invariant. 

The study of the geometry of the bands becomes more useful when the time-reversal 
symmetry is broken. The time-reversal symmetry is either broken by external applied field 



or spontaneously as in [24, |25|. In these cases the Z2 classification is certainly no longer 
applicable. Our natural decoupling of the bands doesn't depend on symmetry of the system, 
so the edge state properties can still be inferred from the study of the geometry of the 
bands. We first consider the case graphene with a uniform exchange field that acts on the 
z component of spin of electron through Zeeman's type of coupling as in 24]. In their work 
the edge states are characterized by the spin Chern number. Though their conclusion is 
qualitatively correct, the physics picture isn't clearly explained and it's somewhat unnatural 
to use spin Chern number when the ground state become ferromagnetic since the spin Chern 
number treats up and down spins equally. In fact the physics picture is quite simple and 
clear in this case. When the Rashba coupling is absent the Hamiltonian can be considered 
as two pieces of Haldane' Hamiltonian with up and down spins. With an increasing positive 
g the relative energy difference of the up spin occupied band and the down spin unoccupied 
band will reduce and the bulk gap vanishes when those two bands meet. If we further 
increase g two bands overlap and the system become ferromagnetic metal. It's the Rashba 
coupling that open a gap at the degenerate points of the two bands which form a surface 
in the three dimensional Brillouin zone and the system become a ferromagnetic insulator. 
From Fig 4 we see the geometry of $jS changes very little with g, it's a reflection that the 
bundle formed by the occupied bands are rather inert in the exchange field as long as the 
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original bulk gap is still open. However, the geometry of the edge state spectrum changes 
a lot, the degenerate cross points are pushed away form time- reversal symmetric point and 
little gaps are opened at the edge. Generally the symmetry-breaking fields may affect the 
edge state spectrum from two different ways, one is to act directly on the edge states, the 
other is by changing the geometry of the bands. The direct action may change the position 
of degenerate points and open gaps at the degenerate point if the fields doesn't commutate 
with the effective edge state Hamiltonian. So even if the indicate gapless edge states, 
gaps can still be opened by direct interaction. In this case, because the geometry structure 
varies very little from the g = case we show the change of the edge state spectrum comes 
solely from the direct interaction. We study the vector spaces spanned by the edge states 
and find there is almost no difference from the g = case. We also compare the edge state 
spectrum and that of a same system except exchange field only present at outmost three 
layers and they are almost identical. When the exchange field large enough the occupied 
bands and the unoccupied bands contact and the geometry structure changes. We see at 
a certain g the two $jS switch partners and the occupied bands naturally decoupled two 
bands each with a Chern number 1 . The system come to a quantum Hall state. In this case 
because of the direct interaction of the large exchange field, part of the edge state spectrum 
have merged into the bulk energy bands. It's no longer easy to compare the geometry of $iS 
and that of the edge state spectrum. However, $jS still predict two branches of edge state 
through the bulk gap at each edge and indicate where the new type of edges state occurs. 



Another important example is the magnetic topological insulators in 25|. Based on 
their first principle calculation result they predict ferromagnetic order appears in the 
Bi2Te3, Bi2Se3, and Sh2Te^ when doped with transition metal elements. Quantized Hall 
conductance can be observed in two dimensional thin film of the ferromagnetic insulator. 
The quantized Hall conductance can be calculated by the Kubo's formula and the change of 
topological structure of the energy bands are discussed in |25]. Certainly we can explain the 
quantized Hall conductance by gapless edge state. We extend their model to a tight binding 
one. When the exchange field is zero in Bi2Te^ film the band structure is just like the 
HgTe quantum well. So we compare the $jS and the edge state spectrum when exchange 
field is nonzero in Fig 6. It's easy to see how change of the geometric structure lead to the 
formation of the quantum Hall type of edge state. The discussion of the change topological 
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structure of the bands in 



25| depend on the fact the Hamiltonian can be decoupled to two 



2x2 blocks and two occupied bands can be easily discriminated. The topological property 
is obtained from the analogy with the Bernevig-Hughes- Zhang model. This property of the 
Hamiltonian can be easily destroyed by e.g Rashba interaction with the substrate and the 
eigenvectors of Hamiltonian won't natually form two bands. In our geometrical approach 
two bands can still be defined and the edge state properties can be analyzed in this way. For 
the same reason the three-dimensional edge state properties of those doped semiconductors 
can't be obtained from the analogy because the Hamiltonian can't be decoupled except at 
kz = 0. So we analysis the edge state property of the three-dimensional material by the 
geometrical property of the bands. Here we use a tight-banding model as in with a 



exchange field in 



25| for doped Bi2Se^ and we consider the edge state in the — plane. 



Without the exchange field the $iS and the edge state spectrums are just like the case 
in k^ — ky plane discussed above except the Dirac cone is anisotropic. When k^ = the 
Hamiltonian has the same form as the two-dimensional case and a large enough exchange 
field can produce quantum Hall type of edge states. When k^ and exchange field is 
zero the $jS and the edge state is gaped. A large enough exchange field can also produce 
the gapless edge state as indicated in Fig 6. However, the farther away from k^ = the 
larger exchange field is needed to produce the gapless edge state. So if a fixed exchange 
field can produce quantum Hall type of edge state at A;^ = when we go far enough from 
kz = the gap will be opened as in Fig 6. At the critical kz = ^kzo the two naturally 
decoupled bands switch partners and become two trivial bands when lA;^! increase further. 
At (0, 0, ±kzo) the occupied bands and the unoccupied bands contact or "kiss" each other 
and form three dimensional anisotropic Dirac points as analogys of the two-dimensional 
Dirac points in graphene. Here through analysis the geometrical structure of the bands we 
predict there is a topological-insulator to topological-semimetal phase transition when bulk 
Bi2Se3 is doped with transition metal elements. From the Fig 6 it can be easily inferred 
there is one Fermi arcs of the edge state as in 261] at each of the opposite surface in the 
topological-semimetal phase. The bulk and surface energy spectrums can be measured by 
ARPES. 



Our calculations with various tight-binding models show that the geometrical structure 
of bands can be a powerful tools to study the the edge state properties of insulators. The 
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algorithm to calculate $iS is very simple. Thus we deem it should be a common procedure 
to calculate the geometry structure of the energy bands in studying the surface structure of 
insulators in the first principle energy band calculation. The geometry structure of the bands 
doesn't depend on the symmetry of the system and provide more edge state information than 
the topological invariants, so our method broadened the scope of the study of the edge state 
properties of the insulators. 
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FIG. 1: Comparison of the <I>jS and the edge state spectrum of some time reversal 
symmetric two dimensional insulators. All pictures in this article are arranged as $jS at left 
and edge state spectrum at right, (a)-(f) show the topological and trivial insulator phase transition 
in graphene when A increases, (g)-(h) show the result in HgTe quantum well. To make the edge 
state more localized we use a parameter M = —0.1. (i)-(j) show the result in a Z2 trivial insulator 
in 
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FIG. 2: Comparison of the <I>jS and the edge state spectrum of some time reversal 
symmetric three dimensional insulators in |l5l| (a)-(b) show the 1; (111) case, there are three 
Dirac points. (c)-(d) show the 1; (111) case, there is one Dirac point, (e)-(f) show the 1; (111) case 
there are no Dirac point, (g)-(h) show the 0; (111) case, there are two Dirac points 
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FIG. 4: Comparison of the and the edge state spectrum of topological crystalline 
insulators (a)-(b) the case discussed in [l5|. (c)-(f) the C4 symmetry is broken by using different 



t2 in X and y directions. In (c)-(f) = 0.7 and = 0-3. In (g)-(h) = 1.5 and = 
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FIG. 5: Comparison of the <I> jS and the edge state spectrum of time-reversal-symmetry- 
broken two dimensional topological insulators. Here parameters of the time-reversal sym- 



metric part of 
is the same in 



lamiltonian in [8j are used with A^, = 0. The Hamiltonian of the exchange field 



24l |. When the parameter g increases the system experience a quantum spin Hall 



state to quantum Hall state phase transistion. 
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FIG. 6: Comparison of the <I>jS and the edge state spectrum of doped Bi2Se^ film and 
three dimensional doped 312833 at different kz (a)-(f) corresponding to the doped 312863 
film or the three dimensional doped Bi28e3 at different kz = 0. With increasing exchange field 
parameter g one of the two bands become trivial, the edge state become quantum Hall type, (g)-(j) 
indicated at A;^ / a large enough exchange field can produce quantum Hall type edge state. For 
a fixed exchange field as in (e) (f) (i) and (j), it's shown in (k) (1) at a certain \kz\ occupied band 
and unoccupied band touch, (m) (n) shows when l^^l increased further the bands become trivial 
and there are on edge states. 
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